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Abstract 

We study a toy model for phantom cosmology recently introduced in the litera- 
ture and consisting of two oscillators, one of which carries negative kinetic energy. 
The results are compared with the exact phase space picture obtained for similar 
dynamical systems describing, respectively, a massive canonical scalar field con- 
formally coupled to the spacetime curvature, and a conformally coupled massive 
phantom. Finally, the dynamical system describing exactly a minimally coupled 
phantom is studied and compared with the toy model. 



1 Introduction 



The idea that the universe underwent an inflationary epoch early in its history [1] has 
spurred interest in scalar field cosmology. A scalar field is the simplest, although not 
the unique, way to fuel inflation. In this context a scalar field, which may come in 
different forms and with a variety of possible self-interaction potentials, is the only or 
the dominant form of matter in the field equations of the gravitational theory. Scalar 
fields have also been considered as candidates for dark matter and, more recently, for 
dark energy or quintessence. The recent discovery that the universe is undergoing an 
accelerated expansion that started not long ago [2] introduces the need to either modify 
gravity or postulate a form of dark energy with exotic properties (negative pressure). 
Scenarios based on modified gravity without extra forms of energy have thus far not 
been very successful. It is easier to describe the observed acceleration of the universe by 
using some form of dark energy, and most models in the literature employ a scalar field 
[3]. 

There is marginal evidence [4] that not only the universe accelerates but perhaps it 
even superaccelerates, i.e., that the Hubble parameter H may increase with time, H > 0, 
which is equivalent to an effective equation of state parameter w = P/p < — 1, where 
p and P are the effective energy density and pressure of the dark energy. If confirmed, 
superacceleration cannot be explained by a canonical, minimally coupled scalar field [5], 
but modifications of the theory are required. One of such modifications, considered by 
many authors [6], consists in a phantom field with negative kinetic energy. Such a field 
is obtained in S-theory, supergravity and in higher derivative theories [7]. The action 
for gravity with a phantom field is 



s (phantom) = J d * x ^j 



(i.i) 



where k = 8irG, R is the Ricci curvature of spacetime [8] and V(4>) is the phantom 
potential. Another possibility to explain superacceleration is to let a scalar field with 
positive kinetic energy couple nonminimally to the Ricci curvature of spacetime, as 
described by the action 



J d^x-J^g 



i^i2-I^v a 0v 6 0-y(0) 



1.2) 



where £ is a dimensionless coupling constant. The explicit nonminimal coupling to the 
curvature introduces extra terms in the equations for the scalar field and the cosmological 
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metric, allowing for the possibility of superacceleration with H > 0. A more general 
coupling between the scalar and gravity, such as the one appearing in scalar-tensor and 
string theories can also generate superacceleration. For example, bouncing universes and 
pole-like inflation, which are not possible in general relativity with a minimally coupled 
canonical scalar, are known in these theories [10]. Quantum effects can also provide a 
mechanism to generate superacceleration from theories that are classically stable (see, 
e.g., Refs. [11, 12]). 

Apart from the need to explain superacceleration, there are other, much more com- 
pelling reasons to consider a nonminimally coupled field. Nonminimal coupling is in- 
troduced by quantum corrections to the action of a classical scalar [13], and is even 
required at the classical level to preserve the Einstein equivalence principle or to avoid 
causal pathologies [14]. 

Of particular interest is the case of conformally coupled scalar field theory with 
£ = 1/6. In fact, this value of the coupling constant is an infrared fixed point of the 
renormalization group [15] and therefore we focus on conformally coupled scalars. 

Once one accepts the idea of a phantom field with negative kinetic energy, it is 
natural to also consider a nonminimally coupled phantom field [16]. While models of 
nonminimally coupled scalars have been studied since the early days of inflationary 
theory, the introduction of phantom fields in cosmology is very recent. The defining 
feature of a phantom, its negative kinetic energy, has been associated with instability 
and runaway solutions. The possibility of negative energies is a well known feature of 
scalar-tensor and generalized gravity and seems to be accepted in this context, at least 
when gravity is classical [10]. Hence, it is not clear that a negative energy is a priori a 
problem. Moreover, when a field is coupled to gravity the gravitational energy plays a 
role in the dynamics and no meaningful prescription is available for the energy density of 
the gravitational field even in general relativity. Issues of principle with negative energy 
are discussed in section 6. Satisfactory answers to the questions raised, however, come 
from the understanding of the dynamics, which constitutes the main part of this paper. 

The problem of the stability of a phantom has motivated the consideration of a toy 
model consisting of two coupled oscillators, one with positive and one with negative def- 
inite kinetic energy, which exhange energy during motion [17]. The canonical oscillator 
with positive kinetic energy is supposed to mimic the gravitational field while the second 
oscillator mimics the phantom field. This dynamical system is also an exact description 
of the dynamics of a massive conformally coupled scalar field in a FLRW universe, using 
suitably rescaled variables and conformal time. Section 2 summarizes the equations and 
motivations for the coupled oscillator model, while section 3 analyzes the dynamics of 
this system. Section 4 considers a similar system of coupled oscillators - both with pos- 
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itive kinetic energy - which describe exactly a conformally coupled phantom in rescaled 
variables and conformal time. Having learned a lesson from section 4, we can approach 
in section 5 the exact description of a minimally coupled phantom, without resorting to 
toy models - the phase space is analyzed in rescaled variables also for this case. Section 6 
contains a discussion and the conclusions. 

2 Coupled oscillators as cosmological models 

In this section we summarize the equations of phantom and scalar field cosmology and 
we show how they can be reduced to a model with two coupled oscillators. We set our 
analysis in the context of a spatially flat Friedmann- Robertson- Walker universe described 
by the line element 

ds 2 = -dt 2 + a 2 {t) (dx 2 + dy 2 + dz 2 ^ (2.1) 

in comoving coordinates (t, x, y, z). In this homogeneous and isotropic universe the scalar 
field (f)(t) only depends on the comoving time and self-interacts through the potential 
V((f>). For simplicity it is assumed that the scalar field is the only form of matter present. 

2.1 Phantom cosmologies 

A minimally coupled phantom field in the metric (2.1) has energy density and pressure 

Pv = -\4> 2 + v X4), (2-2) 

P p = -^ 2 -V(<P), (2.3) 

which correspond to the usual expressions for the energy density and pressure of a 
minimally coupled scalar but with the sign of the kinetic energy inverted. A positive 
potential V(<j>) is necessary if the energy density p p is to be non-negative. Although 
this requirement is sometimes dropped [17] most relativists are reluctant to give up the 
weak energy condition, while all the other energy conditions are progressively being 
abandoned [18]. 

The pressure P p in eq. (2.3) is more negative than for a canonical scalar field, thus 
making it possible to achieve a superquintessence phase with P p < —p p . A minimally 
coupled scalar with canonical kinetic energy density +0 2 /2 in Einstein gravity cannot 
achieve P < —p [5]. 
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The Einstein field equations in the metric (2.1) and with the phantom 
source are 



H 2 = - 



+ 2V((j)) 



H + H 2 = ^ 2 + V{<P) 

■ dV 
+ 3#0- — = 0, 



as material 
(2.4) 
(2.5) 

(2.6) 



where an overdot denotes differentiation with respect to the comoving time t. Only 
two equations in the set (2.4)-(2.6) are independent - when <fi ^ one can derive the 
Klein-Gordon equation (2.6) from the other two, or from the conservation equation 
p p + 3H (P p + pp) = satisfied by the phantom. 

The field equations (2.4)-(2.6) can be derived from the Lagrangian 



Lq = 3aa 2 + kcl 



+ v{4>) 



(2.7) 



or from the Hamiltonian 



Ho = a 3 



nV{<p) 



(2.8) 



A toy model for phantom cosmology consists of two coupled oscillators, one with 
positive-definite and one with negative-definite kinetic energy, described by the La- 
grangian 



x 2 y 2 m 2 x 2 m 2 p 2 

x y xv 

2 2 2 2 y 2 y 



or by the associated Hamiltonian 



x 2 



y 2 



m 



x x 2 + 



m y 2 i A* 2 2 

-y 2/ + Y x y 



(2.9) 



(2.10) 



A similar model was considered in Ref. [17]. However, these authors choose the op- 
posite sign for the term m 2 y 2 /2, and the energy of the phantom turns out to be 

p p = — (0 2 + rri^cf) 2 ^ /2 and negative-definite instead of p p = (— 2 + m 2 (j) 2 ^ /2. While 
in the toy model the negative kinetic and potential energy of the "phantom" oscillator 
y can be offset by the positive interaction energy p 2 x 2 y 2 /2, this is not a possibility for 
the real phantom field. In the exact phantom model, expressions for the interaction 
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energy density between gravitational and phantom fields and for the gravitational field 
energy density are not available, and it is desirable to keep the energy of the phantom 
positive and preserve the weak energy condition. This explains the difference between 
our Lagrangian (2.9) and that of Ref. [17]. 

As remarked in Ref. [17], while the total energy of the system is constant, the energy 
of the phantom oscillator could decrease ad infinitum and the energy of the other oscil- 
lator could increase without bound - the system could not have a stable ground state. 
This argument has been presented before in the context of scalar-tensor gravity, which 
also has negative energies [10]. 

The Euler-Lagrange equations derived from the Lagrangian (2.9) are 

x + (m 2 x + /iV) x = , (2.11) 

y-(m 2 y + f j?x 2 )y = 0. (2.12) 

The physically interesting case corresponds to a massless graviton and m x = 0, while 
the phantom field is very light with a mass ~ 10~ 33 eV. Eq. (2.11) resembles an 
harmonic oscillator equation with a real ^-dependent effective mass, while eq. (2.12) 
exhibits a purely imaginary effective mass i\Jf^1 + y?x 2 . The first property is associated 
with stability and the second one with instability. By changing the sign of the mass term 
for the phantom, as done in Ref. [17], one obtains an effective mass squared (m 2 — fi 2 x 2 ^J 
for the y-oscillator, which is positive and corresponds to stability if fx 2 x 2 < m 2 . Since 
the phantom mass is constrained to be extremely small, one would have to fine-tune the 
coupling parameter fx 2 to achieve stability [17]. 

In the next section we consider the toy model for a massless phantom obtained by 
setting m x = m y = and ji 2 = 1 and described by the dynamical system 

x = - X y 2 , (2.13) 
V = x 2 y, (2.14) 

which is invariant under the transformations (x,y) — > (—x,y) and (x,y) — > (x, —y). 

2.2 Conformally coupled scalar field cosmology 

There are many reasons to believe that a scalar field in a curved space couples nonmin- 
imally to the Ricci curvature: nonminimal coupling is induced by quantum corrections 
to the classical action [13] and is required to renormalize the scalar field theory [19]. In 
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particular, conformal coupling is an infrared fixed point of the renormalization group 
[15] and is required by the Einstein equivalence principle and to avoid acausal propaga- 
tion of 0- waves [14]. Hence, we consider a conformally coupled massive scalar field. The 
equations of motion obtained by varying the action (1.2) with £ = 1/6 and V = m 2 2 /2 
are ([20, 21], [23]-[29]) 

# + 2 # 2 -^V = 0, (2.15) 

2 



+ kH^ - 3H 2 (l - ^0 2 ) + ^- 2 = , (2.16) 

+ 3#0 + -0 + m 2 = , (2.17) 
6 



and they incorporate the Hamiltonian constraint 



H 2 = ^ P . (2.18) 

The effective energy density and pressure of the scalar can be identified in three physi- 
cally inequivalent ways [22] - only one of these guarantees conservation of energy, and 
it yields 

P=Y + ^0 2 + ^0(^0 + 20) , (2.19) 

; 2 2 i 

P = ^L_HL ( p2__ + 20 2 + 200 + (2H + 3H 2 ) 2 ] . (2.20) 

2 2 6" 

The expression (2.19) of p is considerably complicated but the problem can be reduced 
to a system of two oscillators with sharply defined energies in a fictitious Minkowski 
space as follows. By using the rescaled variables [23] -[25] 



' Kin 2 



x = ma , y = y - a (2-21) 

and the conformal time r\ defined by dt = a drj, the equations of motion become 

x" = xy 2 , (2.22) 
y" = ~x 2 y, (2.23) 

where a prime denotes differentiation with respect to rj. Eqs. (2.22) and (2.23) are 
invariant under the transformations (x,y) — > (—x,y) and (x,y) — > (x, —y) and they can 
be derived from the Lagrangian 



The associated Hamiltonian is 

Wi = ^ " ^ + , (2.25) 

and the Hamiltonian constraint is Ji\ = 0. Note that it is incorrect to think of x and y as 
describing gravity and matter, respectively, because y is a mixture of the gravitational 
and scalar field degrees of freedom. However, x = ma is associated with gravity and it 
is its effective kinetic energy that is negative here. 

Formally, the system (2.22) and (2.23) can be obtained from the system (2.13) and 
(2.14) by the exchange of variables (x,y) — > (y,x) and t — > rj. Note however that the 
variables x and y and the time coordinates t and rj used in the two cases have different 
physical meanings. Moreover, for the toy model (2.13) and (2.14) we allow any value of 
x, while in conformal scalar field cosmology the restriction x > applies. 

3 Oscillator model for massless phantom and con- 
formal massive field 

The system of coupled oscillators (2.13) and (2.14) was studied by Castagnino et al. 
[25]. Their results are: 

• The fixed points of the system are (xo,0) and (0,yo) and they are all unstable. 
The system of two oscillators does not have stable positions of equilibrium. 



All the orbits in the phase space (apart, of course, from the fixed points) go to 
infinity as t — > +oo with y(t) — > +oo monotonically while x(t) — > oscillating. 



• There are no cycles and there is no chaos. 

• An asymptotic solution that seems to be generic is 



2 ft 3 \ 

x(f)«-sin(-J , y(t)^t 2 . (3.1) 



The kinetic energy associated with the x-oscillator is 



^(-) = M!~2t 2 cos 2 ^ N ) (3.2) 
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and oscillates with divergent amplitude, while the kinetic energy of the ^/-oscillator, 
K^> = —{i)) 2 /2 ~ — 2t 2 — > — oo. This corresponds to the "instability" described 
by Carroll et al. However, Castagnino et al. argue on the basis of their numerical 
analysis that this solution is an attractor, i.e., it is dynamically stable. 

Technically, while proving that x(t) — > as t — > +00, Castagnino et al. do not prove 
that x(t) oscillates while going to zero. However, this is straightforward to demonstrate. 
Assume that x — > monotonically. Then, if x > the curve representing x(t) has the 
t-axis as horizontal asymptote and x < 0, x > 0. This contradicts eq. (2.13) which 
guarantees that x = —xy 2 < 0. The same reasoning applies if x < goes to zero 
monotonically. Hence, x(t) must oscillate. 

Let us turn now to the case of a massive conformally coupled scalar field. Castagnino 
et al. actually studied the dynamical system (2.13) and (2.14) but their variables are 
not defined as in our eqs. (2.21) - x and y are exchanged in their definition, so their 
system (2.13) and (2.14) is physically equivalent to our system (2.22) and (2.23) and 
their results are derived using conformal time [30]. The proofs in Ref. [25] are then 
given by considering the limit 77 — > +00. This does not necessarily correspond to large 
comoving times, which is what is instead desired. For example, in a de Sitter space the 
scale factor is a = a e Ht = — and t — > +00 corresponds to 77 — > Moreover, if 

big rip solutions of the kind a = (t — t)~ x are present, both comoving and conformal 
time stop at a finite value. However, one can consider the least upper bound 77* on the 
possible values of 77 along the orbits of the solutions and replace the limit rj — > +00 with 
r] — > +77* in the proofs of Ref. [25]. 

When translated in terms of physical variables a and <p and in comoving time t and 
by imposing the restriction x > 0, the interpretation of their results is as follows. 

• The fixed points of the system (x , 0) are Minkowski spaces and they are all un- 
stable. There are no attractor points. 

• All the orbits in the phase space (apart, of course, from the fixed points) go to 
infinity as t — > +00 with y(t) — > +00 monotonically while x(t) — > oscillating. 
These solutions correspond to universes that expand to infinity with scalar field 
that goes to zero oscillating. This result is consistent with the known phenomenol- 
ogy of a nonminimally coupled scalar field: oscillating solutions for £ > (e.g., 
[26]) and non-oscillating ones for £ < (e.g., [27]). 

• There are no cycles and there is no chaos. 
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• The asymptotic solution 

x(t) » 2 - sin ^pj , y(t) » ?? 2 (3.3) 

corresponds to a matter-dominated universe with a(t) = a t 2 ^ 3 . 

An analysis on the lines of that of Ref. [16] leads to the conclusion that the dynamical 
system is not integrable for general values of the parameter \i: however, it is not chaotic. 
(If a quartic self-interaction for the scalar and a cosmological constant are added to the 
scenario, there are special values of the parameters for which the system is integrable 
[28, 24] - however, these cases are fine-tuned and not physically relevant.) 



4 Conformally coupled phantom 

Once one accepts the use of a phantom field, it is natural to consider a nonminimally 
coupled phantom [16]. Here we focus on a massive conformally coupled phantom field 
for which the Klein-Gordon equation is 

+ 3#</> - m 2 - = . (4.1) 

In terms of the rescaled variables x and y of eq. (2.21) and of conformal time r], the 
field equations can be derived from the Lagrangian 

i2 ^ + (^ + I,V (4.2) 

or from the Hamiltonian 

H 2 = ^ + ^- l -xV. (4.3) 

The Lagrangian (4.2) is actually equivalent (apart from an irrelevant overall sign) to 
the Lagrangian L\ of eq. (2.24) in which the conformally coupled scalar is turned into 
a phantom, i.e., in which both gravity and the rescaled phantom have negative kinetic 
energies. 

The dynamical system is [16] 

x" = xy 2 , (4.4) 
y" = x 2 y, (4.5) 
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and the energy constraint is 7^2 = 0. The dynamical system (4.4) and (4.5) is invariant 
under any of the exchanges (x,y) — > (y,x), (x,y) — > (—x,y), or (x,y) — > (x,—y) and 
the phase space is 3-dimensional in these variables. The fixed points are the empty 
Minkowski spaces (x ,0). A series of results analogous to those seen in the previous 
section hold, and they are proven in the same way. 

Theorem all the (nonstationary) orbits in the phase space go to infinity with x — > +oo. 



Theorem there are no closed orbits in the phase space. 



Again, all the solutions that are not stationary represent universes expanding to in- 
finity, all the Minkowski fixed points are unstable, and there is no chaos. Of course, the 
symmetry between x and y leads to features that are different from the previous case: 

Theorem both variables x and y diverge without oscillations as rj — > 77*. 

The proof is as follows. Consider the function I(x,x f ) = x'/x, which is positive along 
orbits different from the fixed points. By differentiating I along the orbits and using 
eq. (4.4) one obtains 

dJ = x 2 y 2 - {x'f ^ ^ 

di] x 2 
The energy constraint H2 = yields 



dl 

drj 




(4.7) 



which is strictly positive along the orbits of non-static solutions. I is positive and strictly 
increasing, hence the limit lim^^ I is either +00 or I > 0. If the limit is +00 then 
^- — > +00 implies that x' — > +00 faster than x, and also x" — > +00. The energy 
constraint implies that 

diverges. Since y 2 — > +00, y cannot oscillate and have zeros. Hence we have one of the 
two possibilities y — > ±00. The second derivative y" = x 2 y has the sign of y and is either 
positive or negative, but it does not change sign. Hence y does not oscillate. 
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If instead the limit of / is finite and there is a horizontal asymptote, then I' = 



(y'/xy -> 0+ Since V > and I' -> 0, it is I" -> 0+. Now, 



(4.9) 



and, by taking the limit of this quantity as rj — > 77* and using the fact that a; — ■> / and 
y'/x — * 0, one obtains 

lim yy' = + . 



(4.10) 



Now we have J" > 0, or yy' > I y^j > 0. Either y > and y' > or y < in 
conjunction with y' < 0. In the first case, in order for y to go to zero, y must decrease 
and it must be y' < 0: we then have a contradiction. The same reasoning applies to the 
case y < with y' < 0, reaching again a contradiction. Hence the limit of / cannot be 
finite. This completes the proof. 

The system is integrable if n — 0, which describes the trivial case of two free non- 
interacting particles. In general, however, the system is not Liouville-integrable and in 
this sense it is dynamically complex [16], but it is not chaotic. 

Finally, we present the two exact solutions 



xi(rj) = yi(v) 



x 2 (v) = -V2(rj) 



±V2 
V~Vo 

TV2 



(4.11) 



(4.12) 



V-Vo 

where 770 is a constant and the sign is chosen so that x > 0. The relation dt = adr) 
is easily integrated to express these solutions in terms of comoving time t and see that 
they represent the expanding or contracting de Sitter spaces 



a(t) 



I m t 

ao e V2 




(4.13) 



This special value of the scalar corresponds to a divergent effective gravitational cou- 



pling K eff = K(l - K(j) 2 /<0) 

universes [29, 33]. 



and to a non-removable curvature singularity in Bianchi I 
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5 Minimally coupled phantom in rescaled variables 



We can now approach exactly the problem of the minimally coupled massive phantom 
without resorting to toy models by using the rescaled variables (2.21) and compare 
the equations of motion with those obtained in the previous sections. In terms of the 
variables x and y, the equations of motion are 



/\2 



x 



X 



+ 3xy 2 , 



x 



2 , 3y 3 + x 2 y , 



while the energy constraint is 



(5.1) 
(5.2) 



(xr+[y'- V —\ =*v, 



X 

and the (rescaled) effective Lagrangian and Hamiltonian (2.7) and (2.8) are 

2" 



1 



X 



x'Y + x 2 y 2 + \y' - 



xy_ 

x , 



X \ x , 



X 



xy 



(5.3) 



(5.4) 



(5.5) 



Eqs. (5.1) and (5.2) are rather complicated in comparison to the cases studied in section 3 
and 4, due to the fact that now the phantom is not conformally coupled. 

The fixed points are again Minkowski spaces (xo, 0). In the new variables we are able 
to derive various results. 



Theorem all the orbits in the phase space (apart from the stationary points) go 
to infinity. 

To prove this statement consider the function F(x,x') = xx' . By differentiating along 
the orbits of the solutions and using eq. (5.1) one obtains F' = 3x 2 y 2 > on all 
the non-static orbits. Since F(rj) is positive and monotonically increasing it is either 
lim^^ F{rj) = C > or lim^^^ t F(i]) = +oo. We show that F cannot tend to a fi- 
nite limit. In fact, if F{rj) — > C, the asymptotic equation xx' ~ C is satisfied and the 
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asymptotic solution is x = \/1Cr\ + D, where D is an integration constant. Then 

C 2 

x" « ^ < , (5.6) 

while eq. (5.1) combined with the energy constraint yields 

x" = -(y'- y —) 2 + 2xy 2 >0, (5.7) 
x \ x J 

in contradiction with eq. (5.6). Hence F = xx' — > +oo, which implies that x — > +oo. 
Theorem there are no closed orbits (cycles) in phase space. 

The proof in Ref. [25] applies also to this case. Hence all the solutions (apart from 
the fixed points) represent universes expanding to infinity, all the Minkowski space fixed 
points are unstable, and there is no chaos. Furthermore, 

Theorem the variable y diverges without oscillating. 

To prove this statement, consider J(x,x') = x'/x, which is well defined for rj — > 77* 
and is positive because F = xx' > and x > 0. By differentiating along the orbits and 
using eq. (5.1) one obtains 

J'(^) = -2^ +3y 2 . (5.8) 
The energy constraint (5.3) then yields 

J'(v) = - 2 (y'-—) 2 + y 2 >o. (5.9) 

x z \ X J 

Hence J is positive and increasing and the limit of J is either finite, lim^^^ J = k > 0, 
or plus infinity. The limit cannot be finite, or else there is a horizontal asymptote of J 

and J' — > 0, which implies that both (y' — ^) — > and y — > 0. But then the energy 
constraint (5.3) implies that 

I X \ If/ ^ 



X / X 2 \ X 



y'-^l +y «_o. (5.10) 
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On the other hand, if J — > k, it is asymptotically x'/x « A; > 0, in contradiction with 
x'/x — > 0. 

Hence, it must be J = x'/x — > +00, which implies that x' — > +00 and, using again 
eq. (5.3), 

y 2 =(-) 2 + - 2 (V - —V > J 2 - +00 (5.11) 

\x J X z \ X J 

and y — > ±00. Moreover, y cannot oscillate or it would go through zeros, which is in 
contradiction with y 2 — > +00. 

Again, one does not expect the system to be integrable, but there is no chaos. 



6 Discussion and conclusions 

We are now ready to summarize and discuss our findings and to compare them with 
previous works. 

The dynamical system (2.13) and (2.14) that constitutes a toy model for a minimally 
coupled massless phantom field and a massless graviton in physical time was studied by 
Castagnino et al. [25]. Their intent was to study the dynamical system describing the 
physically different situation of a conformally coupled massive scalar in rescaled variables 
and conformal time. The latter is obtained from the system (2.13) and (2.14) by the 
exchange (x, y) — > (y, x) and t — > r\. The variables x and y and the conformal time 
employed in Ref. [25] have a physical meaning that is quite different from that of the 
coupled oscillators in the toy model for the phantom. Carroll et al. [17] considered 
the toy model for a phantom but added a negative mass term to the potential in order 
to stabilize the perturbations of the model - the system studied by Castagnino et al. 
corresponds to the special case of that considered by Carroll et al. in which the phantom 
is massless and the perturbations are always unstable. 

The phase space picture of the dynamical system (2.22) and (2.23) describing a 
conformally coupled massive scalar field with canonical kinetic energy density is the 
following [25]. All the fixed points are empty Minkowski spaces and they are unstable. 
There are no cycles and all the orbits of the solutions (except, of course, the fixed points) 
represent universes expanding to infinity, while the scalar field goes to zero oscillating. 
Hence, the matter content (the scalar) tends to infinite dilution in the future while the 
universe expands. This is consistent with the known phenomenology of nonminimally 
coupled scalar fields. 

The phase space of the dynamical system (4.4) and (4.5) describing a conformally 
coupled massive phantom field is qualitatively different: there are again unstable fixed 
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points representing Minkowski spaces and all the other orbits go to infinity, but the 
^/-variable diverges without oscillating. Again, there are no cycles. 

We have studied also the exact phase space of a massive, minimally coupled phantom 
in the rescaled x- and y-variables and using conformal time, without resorting to toy 
models. In these variables, the Lagrangian and Hamiltonian are considerably more 
complicated than in the previous cases, due to the fact that the x- and y- variables are 
suited for conformal coupling. Nevertheless, by using these variables we can compare 
our results for a minimally coupled massive phantom with those for conformally coupled 
canonical scalars and phantoms. The results for this case are as follows. There are again 
Minkowski space fixed points, which are unstable. All the other orbits go to infinity 
with y diverging without oscillations, and there are no cycles. Allowing the phantom to 
be massive does not constrain the orbits to a finite region of the phase space. In fact the 
phantom "falls up" in the potential m 2 (f) 2 /2, which does not have a maximum in which 
can settle. 

There is a difference in the qualitative behaviour of the solutions of the exact equa- 
tions for the phantom in this case and those of the toy model discussed in Refs. [17, 25]. 
While in the exact case studied here x and y diverge without oscillating, the toy model 
has x — > with oscillating behaviour while y — > +oo. The toy model is clearly inade- 
quate and its solutions are wildly different from the solutions of the equations describing 
the actual physical system comprising the gravitational field and the phantom. 

Finally, let us comment on the instability associated with negative energy which 
motivated the introduction of a toy model in the first place [17]. The argument that 
a negative kinetic energy leads to instability is recurrent in Brans-Dicke, scalar-tensor 
and non-linear gravity theories and many authors regard it as a rationale to discard 
the Jordan frame formulation of scalar-tensor theories in favour of their Einstein frame 
counterpart (see Ref. [31] for a review and Ref. [10] for a discussion). This point of view 
is legitimate in flat space physics but is not compelling when gravity is included in the 
picture. First, one cannot write expressions for the energy density of the gravitational 
field (corresponding to the x-oscillator in the toy model) and the interaction energy 
density between the two fields (corresponding to the fj, 2 x 2 y 2 /2 term), and it is not clear 
how such a separation could be achieved. This fact causes ambiguities in the definition 
of energy density for a scalar that couples nonminimally to the curvature through an 
explicit non- minimal coupling term or a Brans-Dicke-like term in the action [22, 10]. 
Indeed, one can associate unambiguously an energy to scalar waves in linearized Brans- 
Dicke gravity, and this energy turns out to be negative - hence one would expect an 
instability. However, the Minkowski background is stable [32] . Second, the discussion of 
runaway solutions associated with a negative energy instability becomes delicate in the 
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cosmological context. A runaway solution representing a universe that expands forever or 
recollapses is undistinguishable from an "ordinary" cosmological solution with the same 
properties, which is regarded as perfectly acceptable. It is true that negative energies 
in scalar-tensor gravity lead to trouble with the formulation of the Cauchy problem or 
with the quantization of the linearized gravitational field [10], but the negative energy 
cosmological solutions seem acceptable at the classical level. 

When the phantom is treated as a quantum field, the instability becomes more 
worrisome because, in general, a quantum field has many more states than a particle and 
there are many more channels for decay - this plays a role in how an instability manifests 
itself. A quantum phantom field can decay into ordinary particles and other phantoms 
and, vice-versa, particles can decay into phantoms and other particles. Since the direct 
coupling of the phantom to ordinary matter must be suppressed to a high degree in 
order not to violate the equivalence principle, there remains the possibility of decays 
of phantoms into gravitons and vice-versa. Since the classical phantom coupled to the 
gravitational field would seem to produce negative energy scalars and positive energy 
gravitons, and gravity couples to the stress-energy tensor of the phantom, the decay 
process would be ultra-fast without a cutoff. Carroll et al. [17] study this possibility 
and conclude that it is not unreasonable to obtain lifetimes larger than the age of the 
universe if a momentum cutoff as low as 10 _3 eV can be motivated. At present, there 
is no justification for such a cutoff, while a natural cutoff would instead be at the 
Planck scale. Hence, although at the classical level the purported instability of the 
classical unperturbed model does not seem to be fatal for phantom cosmology, at the 
quantum level the survival of the universe containing a phantom for 13.7 billion years 
looks surprising. 
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